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Abstract 

The old " glue-and-cut" symmetry of massless propagators, first established in Ref. [lj], leads 
- after reduction to master integrals is performed — to a host of non-trivial relations between 
the latter. The relations constrain the master integrals so tightly that they all can be analytically 
expressed in terms of only few, essentially trivial, watermelon-like integrals. As a consequence we 
arrive at explicit analytical results for all master integrals appearing in the process of reduction 
of massless propagators at three and four loops. The transcendental structure of the results 
suggests a clean explanation of the well-known mystery of the absence of even zetas (C2n) m the 
Adlcr function and other similar functions essentially reducible to massless propagators. Once a 
reduction of massless propagators at five loops is available, our approach should be also applicable 
for explicitly performing the corresponding five-loop master integrals. 
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1. Introduction 

Within perturbation theory quantum-theoretical amplitudes are described by Feynman Inte- 
grals (FI's). The evaluation of the latter has seen quite a lot of progress during last three decades. 
In fact, it has been elevated from a collection of loosely related prescriptions to a solid part of 
mathematical physics as was recently certified by the appearance of Smirnov's bestseller books 
"Evaluating Feynman integrals" and (even!) "Feynman integral calculus" 0, 01- 

A significant number of higher order calculations are performed according to the following 
"standard" scenario. First, the Feynman amplitudes are reduced to a limited set of so-called 
master integrals (Mi's). The particular way of implementing the reduction is not unique and not 
essential for our discussiorQ. At the second and final step the resulting master integrals should be 
computed. 

An important feature of the standard scenario is that the resulting set of master integrals 
should be computed only once and forever due to the well-established^ property of universality: 



1 The so-called Laporta approach [^-[(J seems to be most often utilized but a few other promising methods are 
being now actively developed |7H12|. 

2 At least well-established in practice. See below for an instructive particular example of a class of massless 
propagators and also [13 . for an attempt to formalize the concept of the masters integrals and to prove the 
universality property in general. A related discussion could be found in 
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for every given class of Feynman amplitudes characterized by the number of loops and the pattern 
of external momenta and masses the corresponding set of master integrals is universal in the 
following sense: 

(a) Every (even extremely complicated) amplitude from the class can be expressed in terms of 
one and the same (finite!) set of masters integrals. 

(b) The knowledge of MPs up to some properly fixed order in the e-expansion is enough to 
calculate the finite part of the amplitude. Let us consider an L-loop integral P. The reduction to 
masters leads to an identity of the form: 

P = Y J C l (e = 2-D/2)M l , (1) 

i 

where sum goes over all relevant master integrals and Ci(e) is a rational function of the space 
time dimension D = 4 — 2e and kincmatical parameters likes masses, external momenta, etc. The 
functions could be singular at the point D = 4. The corresponding poles in e are customarily 
referred to as spurious ones. While the coefficients Cj(e) depend, obviously, on the initial integral 
P, the maximal powers, pi, of the spurious poles inside a given Ci(e) depend only on the choice of 
the basis of master integral^. 

Thus, within the standard scenario, to evaluate an L-loop amplitude F one proceeds in three 
main steps: 

(i) Choose a set of master integrals. 

(ii) Reduce every Feynman integral contributing to the amplitude F to form ([I}. 

(Hi) Compute the e — > expansion of each master integral Mi up to (and including) the term of 
order e Pi . 

The steps (i) and (ii) are, in fact, strongly interrelated. In (almost) all approaches to reduction 
one first tries to use the traditional method of Integration By Part (IBP) identified] in order to 
reduce (read simplify) initial integrals as much as possible. The remaining basis set of further 
irreducible (at least in practice) integrals is considered as the set of MPs. As this final set is 
usually rather small it is not of any practical importance whether the corresponding integrals are 
really independent or no1@. 

Once the set of MPs Ms is fixed, then the corresponding powers pi can be easily read off from 
the results of reduction of some test set of initial FPs. Of course if a set of input FPs is too 
limited, it might happen that a in few cases an "experimentally" determined power pi will be 
smaller than its true value. Luckily, the basis set of MPs (together with corresponding maximal 
values of spurious poles in their coefficient functions) is usually determined in early stage after 
calculation of relatively small subset of all FPs to be computed. 

The choice of MPs is not unique. One of the basic criteria is simplicity of the calculation of 
MPs. For example, in view of an analytical evaluation it is natural to seek for MPs with minimal 
number of propagators. On the other hand, for a numerical evaluation it is often advisable to 
consider more complicated but less singular MPs (see, e.g. 0). 

The standard scenario was first developed for massless propagators 21 , 2^, [l| . It is no wonder 
that our understanding of reduction and MPs is most advanced for this case. Indeed, at three-loop 
level there is an explicit algorithm of reduction [l| to MPs (see, Fig. [T]). The existence of such 



3 It was proven in |l9ll that there always exists such a set master integrals that all coefficient functions will be 
regular at e around zero. 

4 In addition to the IBP identities the so-called Lorentz-invariance ones [(J are also often employed in practical 
calculations. In fact, the second set of identities has been proved 1 1 12fl to be a consequence of the first one. 

5 In addition, sometimes there are implicit confirmations of the independence. For instance, if one computes a 
gauge invariant combination of Feynman integrals, then the gauge independence of a coefficient function of a MI 
could be only guaranteed if the latter is independent from all the others, see, e.g. [2o| . 
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an algorithm proves (a)-universality while the (rather tedious) analysis of the structure of the 
algorithm demonstrates that (b)-universality is also valid [lj]. 
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Figure 1: two- and three-loop master p-integrals. e m after a master label stands for the maximal term in £-cxpansion 
of the master integral which one needs to know for evaluation of the contribution of the integral to the final result. 

Let us consider the next loop level in the same class, that is four-loop massless propagators. 
Here the full set of independent Mi's was theoretically constructed in Then a special procedure 
of reduction, based on 1/D expansion of the coefficient functions of Mi's was developed by one of 



the present authors PJ, l23[ with the help of a special parametric representation of FI's, elaborated 
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in [7H9( . The 1 / D method of reduction has been heavily exploited in a series of publications 
in order to compute a number of important physical observables in pQCD. We can not go here 
into the technical details of the four-loop reduction except for the one: it requires huge computer 
storage resources and their effective management. As a consequence its practical implementation 
would hardly be feasible without excellent possibilities for dealing with gigantic data streams 



offered by the computer algebra language FORM [31[ and, especially, its versions ParFORM [32 



3J and TFORM 



Thus, we conclude that the reduction problem for the four-loop massless propagators is solved 
in the practical sense. Analytical evaluation of the corresponding Mi's is the central theme of the 
present work. 

The plan of the paper is as follows. Next two sections provide the reader with general infor- 
mation about the problem. Section 0] explains the essence of our approach in detail on the (now 
easy) example of three-loop master integrals. The really new results are described in section [5] It 
is there we send an expert in multiloop calculations directly. Section [5] discusses perspectives of 
our method as for its extensions to more loops and other kinematical situations. In section [7] we 
demonstrate some peculiar properties of our results which help to solve an old puzzle of absence 
of even zetas from some quantities, like the Adler function, expressible in terms of p-intcgrals. A 
discussion of our results is put in section [5] In the last section |H] we summarize the content of the 
paper and express our gratitude to people and organizations who (which) have been continuously 
supporting us during the painfully long period of preparation of the present publication. 

Our results for all four-loop MPs (together with some auxiliary information) are available (in 
computer-readable form) in 



http : //www-ttp .physik.uni-karlsruhe . de/Progdata/ttplO/ttplO-18. 
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2. Massless Propagators 



Propagators — that is Feynman integrals depending on only one external momentum — ap- 
peared in Quantum Field Theory from its very origin and since then form an important class of 
FFs. Within perturbation theory, every two-point Green function 

G{q)= Jdxe^G(x), G{x) = (0|T \j 2 (x)h(0)] |0), (2) 

with ji and j'2 being in general either elementary fields or (local) composite operators, is expressed 
in terms of propagators. If the momentum transfer q is considered as large with respect to all 
relevant masses, the propagators contributing to G(q) can be effectively considered as massless. 
In what follows we will customarily refer to massless propagator-type FFs as p-integrals. 

p-Integrals appear in many important physical applications. Below we briefly mention some 
most known/ import ant ones (for more details and examples see, e.g. reviews 3^, 33]). 



The total cross-section of e + e~ annihilation into hadrons, the Higgs decay rate into hadrons, 
the semihadronic decay rate of the r lepton and the running of the fine structure coupling arc 
all computable in the high energy limit in terms of p-integrals. This is because these quan- 
tities are either defined in terms of a two-point function <j2j) with properly chosen currents 
ji and j2 or can be reduced to this form via the optical theorem. 

Note, that by high-energy limit we understand not only the case when all masses can be 
neglected but also the possibility to take into account mass effects by exploiting a small 
mass expansion. As a suitable example one could mention the calculation of the power 
suppressed (of order w? q js, m^/s 2 and so on) corrections for the correlators of (axial)vector 
quark currents in higher or ders of pQCD 



Coefficient functions of short distance Operator Product Expansion (OPE) of two composite 
operators can be always expressed in terms of p-integrals with the help of so-called method 
of projectors 43[ 44 1. A good example of an early multiloop OPE calculation is the one 



of the a^. corrections to the Bjorken sum rule for polarized electroproduction and to the 
Gross-Llewellyn Smith sum rule (45| . 

p-Intcgrals are extremely useful in Renormalization Group (RG) calculations within the 
framework of Dimensional Regularization 46-48J and Minimal Subtractions (MS) schemes 
0. 

The naturalness and convenience of the MS-scheme for RG calculations comes from the fol- 
lowing statement jioj |: 

Theorem 1. Any UV counterterm for any FI integral and, consequently, any RG function in 
an arbitrary minimally renormalized model is a polynomial in momenta and masses. 



This observation was effectively employed by A. Vladimirov [5l| to simplify considerably the calcu 



lation of the RG functions. The method was further developed and named Infrared Rearrangement 



(IRR) in [2l|. It essentially amounts to an appropriate transformation of the IR structure of FFs 
by setting zero some external momenta and masses (in some cases after the differentiation is 
performed with respect to the latter). As a result the calculation of UV counterterms is much 
simplified by reducing the problem to evaluating p-integrals. The method of IRR was ultim ately 



refined and freed from unessential complications by inventing a so-called i?*-operation [521 |53|. 
The main use of the R* -operation is in the proof of the following statement {53 1: 

Theorem 2. Any (L+l)-loop UV counterterm for any Feynman integral may be expressed in 
terms of pole and finite parts of some appropriately constructed L-loop p-integrals. 

Theorem 2 is a key tool for multiloop RG calculations as it reduces the general task of evaluation 
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of (L+l)-loop UV counterterms to a well-defined and clearly posed purely mathematical problem: 
the calculation of L-loop p-integrals. In the following we shall refer to the latter as the L-loop 
Problem. 

The one-loop Problem is trivial (see eq. (Q| in the next section). The two- loop Problem was 
solved after inventing and developing the Gegenbauer polynomial technique in x-spacc (GPTX) 



21] . In principle GTPX is applicable to com put e analytically some quite non-trivial three and 
even higher loop p- integral^ (for a review see 57|). However, in practice calculations quickly get 
clumsy, especially for diagrams with numerators. The main breakthrough at the three-loop level 
happened with elaborating the method of integration by parts [i^, El of dimensionally regularized 
integrals. All (about a dozen) topologically different families of three-loop p-integrals were neatly 
analyzed in [l| and a explicit calculational algorithm was suggested for every case. As a result the 
algorithm of integration by parts for three-loop p-integrals was established. Later the algorithm 
was implemented (and named MINCER) within the computer algebra languages SCHOONSCHIP 
58| and FORM ^ (see Refs ^ and respectively). The most recent FORM version of 
MINCER is freely available from http : //www . nikhef . nl/~f orm. 

During last two decades MINCER has been used intensively to perform a number of impor- 
tant calculations of higher order radiative corrections in various field theories. As a couple of 
outstanding examples, characterizing the issue, we mention the analytical evaluation of the 0(aj?) 
correction to the ratio R in massless QCD 61, 62| and recent analytical calculations of three- loop 
deep-inelastic structure functional [63l - J65l |. 

Note, that every L-loop Problem is naturally decomposed in two: (A) reduction of a generic 
L-loop p-integrals to masters and (B) evaluation of the latter. As A-problem has already been 
discussed, we proceed now to B-problcm. For L equal to 1 or 2 problem B degenerates to a trivial 
one due to the fact that all masters, being primitive ones, are easily evaluated in terms of T- 
functions. At three-loop level there exist only two non-triviaJf| master integrals whose evaluation 
was rather simple with the help of GPTX and, in fact, was performed well before the algorithm of 
reduction of three-loop p-integrals was discovered. Thus, in three-loops A and B problems could 
be considered as two separate ones. 

The situation is different in four loops. In this case there exist 14[ twenty eight master integrals 
pictured on Fig. [2] and only 15 of them (all after M43) are simple. We call a four-loop p-intcgral 
simple if it is either primitive or reducible to the so-called generalized two-loop F-diagram with 
insertions, F{n\ + a\e, . . . ) pictured on Fig. [3] The corresponding F-integral has been intensively 
studied since work [6(1 ] and is now in some sense analytically known 67 1. The remaining 13 masters 
(from Mqi till M43, reading the table from left to right and the from top to bottom) happen to 
be quite difficult to deal with even numerically, not speaking about analytical evaluation. 

The aim of the present work is to demonstrate that there exists a remarkable bootstrap-like 
connection between parts A and B for the L-loop Problem irrespectively the specific value of L. 
The connection is powerful enough to result in an explicit solution of problem B for L equal to 
three and four (which we will demonstrate explicitly) and, in all probability, for L equal to five 
(we will provide the reader with a strong argument for it). 

The only prerequisite for our considerations is the solution of problem A for the corresponding 
number of loops. 



6 The GPTX is also ideally suited for high-precision numerical calculations of finite p-integrals (with simple or, 
better, without numerators) with really many loops. See l5445dl for a number of spectacular examples in four, five, 
six and even seven loops. 

7 These calculations, in fact, have required development and application of a number of additional technical tools 
(including highly advanced version of integration by parts algorithm) than just the use of MINCER and the method 
of projectors; please consult the original works for further details. 

8 More precisely, we mean non-primitive p-integrals, see definitions below in section l3l 
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Figure 2: all master p-integrals for the four-loop Problem. In Mij the digit i stands for the number of (internal) 
lines in the integral minus five and j numerates different integrals with the same value of i. The integrals are ordered 
(if read from left to right and then from top to bottom) according to their complexity. e m after Mij stands for the 
maximal term in e-expansion of Mij which one needs to know for evaluation of the contribution of the integral to 
the final result for a four-loop p-integral after reduction is done. In other words, m stands for the maximal power 
of a spurious pole l/e m which could appear in front of Mij in the process of reduction to masters. 
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Figure 3: the generalized two-loop p-intcgral; indexes besides lines show the powers of corresponding masslcss 
propagators. n$ and are assumed to be integers. 



Wc will describe how the good old " glue-and-cut" symmetry of masslcss propagators [21[ leads 
— after the reduction to master integrals is performed — to non-trivial relations between the latter. 
The relations constraint the masters integrals so tightly that they can all be analytically expressed 
in terms of only few, essentially trivial, watermelon- like integrals (see diagrams M31, M01, M12, M\\ 
and M23 on Fig[2]) . This provide us with explicit analytical results for all master integrals appearing 
in the process of reduction of masslcss propagators at three and four loops. By an analytical result 
we mean, of course, not an analytical expression for a master integral taken at a generic value 
of the space-time dimension D (which is usually not possible except for the simplest cases) , but 
rather the one for proper number of terms in its Laurent expansion in D around the physical value 
D = 4 as it was discussed in detail above in section [T] 

Note, that for our aims it is completely irrelevant how exactly the part A (reduction to mas- 
ters) is performed/implemented. In fact, we only need the reduction for relatively simple cases 
of p-integrals: namely, no squared propagators and relatively low powers of scalar products in 
numerators. In particular, no knowledge of (admittedly rather complicated) reduction techniques 
based on the asymptotic 1 /D expansion is necessary For understanding of all considerations of 
the paper it is enough to assume that the reduction (problem A) is done with some implementation 
of the Laporta algorithm. 

3. Recursively one-loop integrals 

Without loss of generality we will consider the scalar p-integrals defined in Euclidean space- 
time. Let F(q,e) be a dimensionally regulated scalar L-loop p-integral depending on external 
momentum q and the space-time dimension D = 4 — 2 e. Its dependence on q can be written as 

F(q,e) = f(e)(q 2 r/ 2 -^ (3) 

where uj is the canonical mass dimension of F(q, 0) and f(e) is a meromorphic function of e. 

The complexity of computing of the function of f(e) depends on the loop number L. At one 
loop level the result for the generic integral 



/ 



{P a ){q-l)W 
is known since long (see, e.g. 0) and reads@ 



( 9 2) 2 -«— ?G{a,P) (4) 



r(a + /3-2 + e)r(2 -a -e)r(2 -/?- £ ) 
G{ °' P)= r(a)r(/3) r(4-a-0-2e) ' (5) 



9 We provide every loop integration dPl with an extra normalization factor l/7r D / 2 and write • 
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Here "generic" means that the powers a and (3 could be not only integers but functions of e. The 
most useful in applications case is 



a = m + ae, {3 — n + be (6) 

with n, m being arbitrary integers, a,b nonnegative ones. Note that negative values of a and/or b 
might lead to e independent singular factor(s) like r(0) within the corresponding G-function. On 
formal grounds G(a, /3) is not defined in this situation 10 !. The reduction formula for G-functions 

G(a, ft) = W - W - J — G(a, ft 1) (7) 

as well as the expansion 

G(l + ae,l + 6e) = G ° (e) f 1 + (a + b)e + (a + 6)(a + 6 + 2)e 2 + ■■■), (8) 
e(l + a + b) V / 

G (e) = eG(l,l) = 1 + e (2 - lE ) + . . . (9) 

allows for a convenient evaluation of G(n + ae, m + fee) without any reference to the awkward 
formula (|5|). In fact, the well-known freedom in the definition of the dimensional regularizatioil^l 
allows to tune the function Go(e) at will (provided Go(0) = 1). The most natural choice 

G(l,l) = i (10) 

or, equivalently, 

G (e) = 1 (11) 



fixes the so-called G-scheme [21( and will be adopted here. Note that the G-scheme is not only 
extremely convenient from purely calculational point of view; it is also "natural" in the realm of 
masslcss propagators. There is evidence that results expressed in the G-scheme usually tend to 
display a better pattern of "apparent" convergence in comparison to the MS scheme. 

In view of eqs. © and ^) any recursively one-loop p-integral can be easily performed analyti- 



cally [21(. We will denote such integrals primitive ones. For example, the two- loop Mi's T\ and Ti 
(see Fig. 1) are both primitive ones, their e-expansions (with accuracy necessary for the-two-loop 
calculation) can be easily computed via G-functions: 

Ti = -J--5-^ +£2 p^ + ^W3 ) , (12) 
4e 8 16 V 32 2 / 

T 2 = \ + 0{e 2 ), (13) 
e^ 

with C„ = X);>i 7^"- Here and almost everywhere below we set q 2 = 1. 

For future reference we provide below expressions in terms of G-functions for the four watermelon- 
like primitive three-loop master integrals which serve as building blocks for all other (three-loop) 
masters (see section^. To make the formulas shorter we always use the G-scheme defining relation 
(fTUl) and write everywhere 1/e instead of the G(l, 1): 

P 1 = ^G(2e,l), P 2 = ±G(e,l), 

P 3 = ^G(e,e), P * = 7E- ( 14 ) 



10 See, however, 16811 for a significantly deeper discussion of such cases. 

11 The freedom amounts to the multiplication of every L-loop integral by a factor n(e) L , with n(e) = 1 +£>(£) being 
a regular (at least in a vicinity of the point e = 0) function of e |6fll . Thus, the formulas (1 1 1 and (TTTJ below should 
be understood in the sense that n(e) is chosen as follows n(e) = l/(e G(l, 1)) = T(2 - 2e)/(r(l + e)r(l - e) 2 ). 
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4. Three-loop integrals 



In this section we discuss the main idea of our method on a first non-trivial example of three- 
loop massless propagator-like integrals. 

4-.1. Three-loop finite p-integrals and glueing 

It is easier to explain the glue-and-cut symmetry on a real-life example. Almost exactly thirty 
years ago one of the present authors was strongly puzzled by the following facts (resulting from 
first calculations made with the help of just discovered technique of Gegenbauer polynomials in 
the position space 2l|): 



L = 

JVi = 



(q 2 )- 1 ^ 20 Cs 



0{s), 
O(e), 



N = (qT'~ 6£ 20 Cs 
N 2 = (q 2 )- 1 - 3 " 20 Cs 



O(e), 
0(e), 



(15) 



where L ,N ,Ni and N 2 are scalar three- loop p-integrals (see Fig. 2]). 




.?=() 

9^1 




_e=0_ 





20 Cs 



L N JVi N 2 

Figure 4: Four finite three- loop p-intcgrals displaying a remarkable feature of being equal at e = and q 2 = 1. 



Indeed, a short look on cqs. (|15[) immediately leads to an obvious question: why on the Earth 
four quite different (but all finite) p-i nteg rals have identical values at D = 4 (if one set q 2 = 1)? 
Incidentally, by that time a pioneering^ calculation of the four-loop /3-function in the </> 4 -model 
[7oj | had been just finished. One of its results was the UV divergence of the following four-loop 
vertex-type integral 



UV 




5Cs 

£ 



(16) 



The suspicious appearance of one and the same irrational constant Cs with very simple coeffi- 
cients in eqs. (|15[) and (|16[) was suggesting some mysterious connection between three-loop finite 
p-intcgrals Lq,Nq, Ni, N 2 and the divergent part of the four-loop vertex-type integral E4. In 
addition, a closer inspection of all five diagrams revealed that the four propagators-type diagram 
could be formally produced from the vertex graph in two steps (see Fig [SJ : 

(i) Delete all four external lines from the vertex diagram (transforming it, thus, to a vacuum one). 

(ii) Cut in the resulting vacuum diagram either a line (there exist only two non-equivalent choices, 
leading to iVi and N 2 ) or delete the central verteJ^I (again, one could do it in two ways, leading 
to A^o and Lq). 



12 To our knowledge it was the first full calculation of the /3-function in a four-dimensional model in four loops. 

13 The operation of deleting of a vertex means that one first transforms the vertex into two new ones by introducing 
a fictitious line (with the unit propagator) and then cutting the new line. Note that deleting a three-linear vertex 
does not produce any new diagrams in addition to those coming from cutting the corresponding three incident 
lines. In general, one can cut a four-linear vertex by a three non-equivalent ways a shown by Fig[5jd,e,f). Due to 
the high symmetry of the envelope E4 diagram the possibility (d) and (e) lead to one and the same result. 
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(a) 



(b) 







(c) 



(d) 



(e) 



(f) 



Figure 5: Two ways of cutting a line in the graph E4 (a,b). The generic four-linear vertex (c) and three ways 
transforming it into a pair of the three-linear ones connected by an auxiliary propagator (d,e,f). 



The puzzle was finally understood after the geometrical construction was provided with an- 
alytical content. As a result the Glue-and-Cut (GaC) symmetry of masslcss propagators was 
established. Below we prove a theorem [l[ which solves the puzzle. 

Let (r)(e) be a dimensionally regulated massless scalar L + l loop vacuum Feynman amplitude 
without any subdivergences and with the superficial divergence index wr = at four-dimensions. 
Surely, every expert would cry at this point that such an object is identical zero due to absence 
of any intrinsic scale which is true, beyond any doubtJ^l. Please, be patient! In fact, by a 
Feynman amplitude we understand a formal triplet consisting of the corresponding Feynman graph 
r, properly constructed Feynman integrand and, at last, a function of kinematical parameters 
(external momenta and masses) resulting after evaluation the integrand. It means that even if the 
function vanishes for some particular choice of the kinematical parameters it may become nonzero 
after some modification of the latter. 

Without essential loss of generality we assume that the graph V contains only triple vertexes. 
Consider an arbitrary line of £ of T with 



Pe(q) 



being the corresponding propagator of T and V(q) being some polynomial in the line momentum q. 
We allow the integral (r) to contain non-trivial numerator; in the case of the line I being a fictitious 
one with the unit propagator, one can always redefine the propagator as follows: 



V\q) 



with V'{q) = q 2 



Let (r)(mo,e) be the Feynman amplitude obtained from (r)(e) by introducing an auxiliary non- 



14 By the way, a real mathematical justification why one could self- consistently set zero such massless vacuum 
integrals within dimensional regularization was to best of our knowledge performed only in l7lll . It requires first 
of all a self-consistent and mathematically solid definitions of the very dimensional regularization which, in turn, 
demands a heavy use of various parametric representations. 
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zero mass mo into the selected propagator, as an infrared regulator, that is, 

a 2 Via ) 

Pi(q) -> Pe{q) = ; q :~ J - 

q 2 + m,Q q 1 + 

Assuming that the loop momenta in (r) are chosen in such a way that the momentum q is the 
loop one, we can formally present the integral (T)(mo,e) as a convolution 

<r>(m , £ )= f{rV)(q,e) J^-*, (17) 



where (T\£)(q,e) is the L-loop p-integral obtained from (r)(mo,e), first, by "freezing" integration 
over q and, second, by multiplying the result by (q 2 + mfy. 

Theorem 3. Under the above listed conditions the following statements are true 

(a) The vacuum integral (r)(m-o, e) is IR finite and its UV divergence is a simple pole, that is 

<r)(m ,e) = — ^— + 0(e°), (18) 
[L + l)e 

with C being a constant; 

(b) for every choice of the line I the p-integral (r\^)(g, e) is finite and its value at D = 4 meets 
the condition 

]im(T\£)( q ,e) = (19) 
£->0 q z 

Proof. 

(a) Due to the assumed absence of any subdivergences, the insertion of a mass docs not influence 
the UV divergence, but removes the only possibility for the IR one (related with integration 
over the region of all loop momenta being small) . On dimensional grounds we have 

{r)(m , e ) = K)-^+ 1 ) £ /(e) J (20) 

with /(e) depending only on e. Without any UV subdivergences, the (minimal) UV coun- 
terterm corresponding to the integral (G)(mo,e) as a whole reduces to its pole part. If /(e) 
would contain a non-simple pole it would lead to appearance a non-polynomial dependence 
on the mass mo of the countertcrm in the direct violation of Theorem 1. Thus, 

C= (L + l) lime /(e). (21) 

(b) The (scalar) p-integral (T\£) has a homogeneous dependence on q 2 , namely 

(T\e}( q ,e)^( q 2 )- 1 - L£ g(s) (22) 

with 17(e) depending only on e. As a consequence, the integral over q in (fl~7|) can be easily 
performed with the help of a textbook formula (see, e.g. [i|): 



(m 2 +£ 2 )(e 2 ) 1 + n£ 

(m , )D /«« r((n + i) e) r(i („ + 1 ) e) = 1 

1 (2 — e) [n+ l)e 

with the result 

(T)(mo,e) = -0^- £ +O(e°). (23) 
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A comparison of eqs. (|20I21[) and ((23]) directly leads to eq. 



lim q(e) = C 

which is equivalent to eq. ([lU|) . □. 

The GaC symmetry, proven in Theorem 3, clearly explains the origin of relations displayed 
in Fig. |U Still, considered by itself, it is not especially useful as it does not provide us with the 
value of the constant C. 

4-2. Three-loop master integrals from glueing 

The situation is radically changed if one utilizes the GaC symmetry together with the reduction 
to masterJ^l. Indeed, let us forget for the moment about eqs. (fl~5j) and use only the GaC symmetry 
for the p-integrals shown on Fig. 2] . This leads to four equations, namely: 

N = L + 0(e), N = N 1 + O(e), N = N 2 + O(e), N = O(e ). (24) 

On the other hand the reduction of three reducible p-integrals in eq. to masters gives: 

_ 3 (3D -10) (D - 3) 4(D-3) 2 32 (2D - 7) (D - 3) 2 

L ° (D~ 4^ (D-W 4+ 4p Pl (25) 

12(3D - 8) (3D - 10) (D-3) 4(9D 2 - 65D + 118) (3D - 8) (2D - 5) 

(D - 4)3 2+ (£> - 4) 4 3 ' 

(3D - 10) (D - 3) 8(2D-7) (D - 3) 2 

^ = (D~ 4p il+ (b—^ Pl (26) 

2 (3D - 8) (3D - 10) (D — 3) 4(3D - 8) (2D 5) 

■ "2 + ~ — "a 



(D-4) 3 z (D-4) 5 

16(2 D-7) (D-3) 2 (3D -8) (3D - 10 ) (D - 3) 
(D - 4)3 1+ (D - 4) 3 

10(3D - 8) (3D - 10) (2D - 5) (2D - 7) 



^ /r> A\3 Pl + (n-A\3 P2 V 1 ' 



(D - Af 



D3 



Now, it is well-known fact that a maximal order of the pole in s of a (dimcnsionally regulated) 
L-loop p-integral can not exceed^ L. Thus, we can parametrize the coefficients of the three-loop 
master p-integrals as follow^: 



No = J2 N ^ £l + °( £ )> X i = J2 L M £l + 

i— — 3 

3 3 

Pi = E p i» £l + °( £4 )' P2 = E p ^ £l + 



i=— 3 i— — 3 

-1 



p 3 = E ^W + 0(e 5 ), D 4 = %^, P 4l - 3 = 1. (28) 



i=-3 



Note that the higher term in e, not shown explicitly on (|28p. at any case can not, obviously, be 
constrained by eqs. (|28[) . In addition, they could contribute only to terms of order e or higher 



15 Unfortunately, this step was overlooked thirty years ago, presumably, because the problem of evaluation of 
three-loop p-masters had been already solved before the idea of glueing appeared. 



16 In l53l the statement was proved for an arbitrary euclidean Feynman integral. 

17 Below the explicit result for the simplest master integral, P4 is taken as granted; this fixes the global normal- 
ization of all remaining integrals. 
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to the value of an arbitrary three-loop p-integral (this statement follows from the knowledge of 
maximal power of spurious pole in e which might appear in front of a master integral in the process 
of reduction, see Fig. [T]). 

After substitution of eqs. (|2"5)l , eqs. (f2"4"|) produce some non-trivial constraints on coefficients 
of the e-cxpansion of our master p- integrals. To be specific, let consider first eq. from (12"4"1) . Its 
expanded in e form reads 

lim£™(L - iV ) = (29) 

for n > 0. Note that ([29| is met identically if n > 7 (because of the fact that the maximal allowed 
poles in e which could appear in eqs. (I25I26I27P and eqs (f2"5|). are 4 and 3 respectively). For 
n = 7, 6, 5 and 4 the resulting equations are, 



18 



= P 3 ,-3, (30) 

= 6P 3> _2 + 12P2,-3-4Pl J _ 3j (31) 

O T 

= -59P3,_2 + 6P 3) _i-78P2_3 + 12P2,-2 + 32P 1 _3-4Pi,_2 + ^ ;:i + l, (32) 
= 239 J P 3 ,_2-59P 3 ,-i + 6P3,o + 162P2,-3-78P 2 ,_ 2 + 12P 2 ,_ 1 (33) 

-80 P Xi _ 3 + 32 Pi _ 2 -4P W - ^± + ^ - 4. (34) 

Already now we can see that eqs. (|3"2"f and (|34p express two coefficients, £1,-3 and Li,-2) of a 
non-primitive integral through the coefficients of primitive ones, namely, Pi, P 2 and P3. Indeed, 
a solution of eqs. (|30H34[) is 



p3,-3 


= 0, 








(35) 


Pl,-3 


_ 3P 3 ,_2 

2 +m - 3 > 








(36) 


L1.-3 


22P 3 _ 2 

= 3 4P3 — 


12P2,-3 8P2 r 2 +g ^ , 






(37) 


L\,-2 


128P 3 _ 2 58P 3 
= tut _| i 

3 3 


- - 4P 3 , - 8P 2 ,_ 3 + 12P 2 ,- 2 - 





2 

~ 3 


(38) 



Proceeding in the same vein we arrive eventually to a linear system of 31 equations (not neces- 
sarily independent) for 32 coefficients No t i , 1 Pl.ji j P2J2 > P3J3 ■ O ne can solve the system by 
expressing the coefficients from more complicated masters through those from less complicated 
ones. A convenient ordering is given by two ruleJ^I 



• N 0tlQ >- L Ml y P hjl y P 2j2 y P 3j3 y P 4 



.1 1 ■ 



• For two coefficients of a master integral the more complicated one is that with larger value 
of the second index. 

The system is easily solved with the result: coefficients are expressed in terms of only eight 
coefficients of three primitive integrals, namely, 

P3.-1,, P3,0, P3.I, P3.21 P3,3, P3,4j P2,3j Pi, -3- (39) 

These eight coefficients are trivially determined from eqs. (|14[) and, finally, we arrive at the fol- 
lowing results for the three-loop master integrals: 



18 Note, that for brevity in writing eqs. H3H34I I we have used eq. 1301 to discard the terms proportional to (zero) 
coefficient P3 : —3. 

19 The order between primitive integrals is rather arbitrary, except for the natural choice to use P4 as the easiest 
one. 
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= 20Cs+ 0(e), 



(40) 



111 7 14( 3 / 67 14C 3 ^ 

+ £2 ("IT + ^ + 7C4 + 126 Cs ) + ° (e3) ' (41) 

+ £ 3 (-^ + ^ + 32 C4 + 64 C5 ) + ^), (42) 

J_ A_^Z + f_l^ + i^U 2 / 891 15C 3 9C 4 \ 
4e 2 8e 16 V 32 2 / V 64 4 4 / 



5265 81 Ca 
1 5i _ 


45 Ci 


21 Cs 


128 8 


h 8 ~ 


2 



1 


1 


35 




991 


e+e 2 ( 


26207 


11 Ca 


36 


- + 

e 


216 


+ 


1296 




7776 


18 



(44) 



670951 


385 Ca 


11 & 


46656 


108 


12 



15 




16852031 10901 Cs 385 C 4 13 Cs 



279936 648 72 2 



) 



(45) 



<xx> 



(46) 



where we have boxed the eight input coefficients. The comparison with eqs. (|14I15[) and (|25l) 
demonstrates the all unboxed coefficients have been correctly determined through the gluing pro- 
cedure. 

A remarkable feature of the above discussed, glue-and-cut based, determination of the three- 
loop masters is that the both non-primitive (read non-trivial) master integrals No and L\ have 
been expressed through essentially trivial (read primitive) FI's. Even more, as many as eleven 
coefficients of the primitive Mi's P\ and P2 (see eqs. (|42l43p ^ have also been fixed through only 
eight coefficients listed in ([39]) . Thus, we see that integration by parts identities together with the 
glue-and-cut symmetry severely constrain the values of master integrals. 

5. Four-loop Integrals 

5.1. Four-loop master integrals from glueing 

Following the same procedure in the case of four-loop propagator massless integrals, one should 
consider all possible cuttings of a set of five-loop vacuum massless diagrams with integrand of mass 
dimension twenty and without subdivergences (or, equivalcntly, superficially and logarithmically 
divergent). 

Again, as in three-loop case, GaC relations provide us with enough information to express 
all the necessary coefficients of the e-cxpansions of all MPs through some trivial integrals. The 
number of the input five-loop tadpoles and the resulting relations (around a hundred and a 
thousand respectively) are too large to be presented here, so in the equations to follow we give 
only the final results. 






(47) 




(48) 



M 62 , e 1 
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A/51, e 1 

. £ (-195C 5 + 153C! + ^-51( 



C4 - ^) + 0(e% 




^£ -80C 5 -22C| + 50C 6 



8OC5 + 88 c 3 2 ~ 200C 6 - 66 C 3 C 4 + 




20C5 



8OC5 + 8C3 +50Ce 



M 42 , e 1 

e (8OC5 - 32 Cl - 200C 6 + 24 C 3 C 4 + 520 (7) + 0(e 2 ), 




441 C? 



0(e 



^44. 




= 36 Cf + e (108 C3 C 4 - 378 Ct) + 0(e 2 ), 
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1 7 1 / 17 25Cs\ 377 25C 3 25( 4 



12e 4 4e 3 12e 2 s V 12 6 7 12 



M 34 , e 3 

3401 463C 3 75C 4 465C 5 



e 



e 



12 6 4 2 
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12 6 4 2 

2 / 24497 3031C 3 463( 4 1395C 5 1247(1 3425Ce 



12 6 4 2 6 6 

3 / 158273 19663C 3 3031C 4 6807Cs 
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3(3 

2 
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437C 5 



87C| 



- 90Ce + 
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20C5 



80C 5 + 68C 3 2 + 50C 6 



M 52 , e 1 
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■^+45C 5 -17CI-^ 



'43: 
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e [ -195C5 + 153 C 3 2 + ^ - 51 C 3 C4 ^ ^ ) + 0(e% 



(59) 






1 1 


-H 






67 


) 3e 4 


3e 3 1 3e 2 




-5- 1 
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M 32 , £ 3 
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M 3 i, e 3 

Thus, we observe that at the four-loop level the GaC method works as good as the three-loop 
ones: all required terms of the e-expansion of every four-loop MI have been expressed in terms 
of only twelve coefficients (boxed in eqs. (|47ti74[) 1 

M 23 ,4, M U ,3, M U ,4, Mii, B , Mi 2>4 , Mi2, 5 , M 1,2, M 1,3, Af 01 ,4, M 1, 5 , Moi,6, Mai -4, (75) 



of primitive watermelon-like massless propagator integrals. 
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An inspection of the above results for Mi's reveals a few remarkable features. 

1. In agreement with common expectations (based on the known solutions of the two- and 
three-loop B-problem) the transcendental terms up to (and including) weight 7 appear in 
cqs. (|47H73[) . That is all results depend on only five irrational constants: £ 3 , £4, £5, and 
Cr. 

2. For a given MI Mj the term e Pi (that is one with maximal power in e) aZwayj^l includes £7. 

3. A term proportional to e Pi ~ 3 could contain Q n with n not exceeding 7 — j; if 7 — j < 3 then 
the term is free from irrational numbers. 

4. There is another restriction on the singular part of any MPs (in fact, it is valid for an 
arbitrary p-intergral). It states that the term e~ n (with n = 1,2, 3,4) may not contain zetas 
with the transccndcntally weight exceeding (7 — 2n). This property explains a very peculiar 
feature of MPs M§2 and Ms 3 : the absence of m the corresponding 0(e Pi ^ 1 ) terms. 

5. The only two finite MPs, namely, M44 and M45 contain only terms of one and same weight 
in every (available) coefficient of their e-expansions. 

6. The same property of "transcendental homogeneity" is true for the MI M52 (which is up to 
a factor of 1/e is the three-loop finite MI Nq) if one divides an extra factor (1 — 2e) 2 out of 



it. (See in this connection work 72j, where some general arguments were given in favour of 



the hypothesis that the property is valid in all orders in e.) 

We want to stress that any statement on the structure of £'s appearing in an integral does 
depend on the global normalization which is rather arbitrary. Our normalization condition is a 
natural one but, certainly, not unique. If we would choose 

m 3 i = 4(1+ J2 a * £i ) ( 76 ) 

1<?<7 

then all MPs would depend on aj, and all statements just discussed above could be, obviously, 
made invalid if a coefficient at were allowed to contain Q (for i > 1) and je for i = 1. 

On the other hand, if we restrict ourselves to a natural choice of the normalization of MI M31 
such as 

M 3 i = ^(l+ (77) 

3<i<7 

with bi being rational numbers, then the properties 1-6 would in general stay untouched. 

5. 2. Tests of the results 

In this subsection we discuss various checks which we have made to test our results expressed 
in eqs. (|47H74j) . The set of 28 master integrals is naturally divided in three subsets: primitive 
(M 23 , M 2i , M 25 , M u , Mia , M 13 , M u , M i , M 31 ), simple (M 32 , M 33 , M 21 , M 22 , M 26 , M 27 ) and, 
finally, complicated ones (M 6i , M 62 , M 63 , M51 , M41 , M 42 , M 44 , M 45 , M 34 , M 35 , M 36 , M 52 , M 43 ). 
We will consider these subsets separately. 

5.2.1. Primitive integrals 

A primitive FI is by definition expressible in terms of the r-function. A straightforward use of 
formulas of section [3] gives: 



20 MI M31 is the only exception from this rule since its sub-leading in e terms are fixed to be zero essentially by 
hands, that is by choosing the G-schcmc. 
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M 23 = ^G(e,l), M 24 = ^ G(2e, 1), M 25 = 1 G(3e, 1), 

Mn = -^G(e,e), M X2 = -i G(e, 1) G(e, 1), M 13 = G(e, 1) G(-l + 3e, 1), 

M 14 = ^G(e,2e), M i = ^ G(e, e) G(-2 + 3e, 1), M 31 = 1. (78) 

After the expansion in e eqs. ((78)) produce altogether 89 coefficients. As was discussed in 
subsection 15. II as many as twelve coefficients listed in eq. (|75|) have been used in the process of the 
solution of the system of GaC equations, while the remaining 77 coefficients have been predicted 
from the equation and listed unboxed in eqs. (|66H73I) . 

The reader is advised to check that all these 77 coefficients are in full agreement to eqs. (|66t473p ■ 

5.2.2. Simple integrals 

These all could be expressed in terms of G-functions and the generalized two-loop diagram: 



M 32 = -F(l,l,l,l,e), (79) 

M 33 - ^F(l,l,l,l,2e) (80) 

M 2 i = ^F(l,l,l,e,e), (81) 

M 22 = 4rF(l,e,l,£,l), (82) 

M 26 = Jg(3£,1)F(1,1,1,1, £ ), (83) 

M 27 = 1 0(6,1)^(1, 1,1,1,2s-!). (84) 



To be specific, let us consider the direct evaluation of M 33 in some details. First, we define a 
related FI M 33 pictured in eq. (|85|) below. A simple reduction of FI M' zz to MPs gives the (exact) 
equation 



<2> 



4(7/? - 26)(5 J D - 14) (5D - 16)(D -3) -2(2£> - 7) 

Mi4 H ^7-^ — M 33 . 



9(3D-10)(£)-4) £ 



3(£>-4) 



^33 



On the other hand, 

M^ 3 = lF(l,l,l,l,l + 2e) 

and (see [z3 as well as HSIfi^l) 
1 



F(l, 1, 1,1,1 + 2s) = 6C 3 + 9C 4 e + 192C 5 e 2 + (465C 6 - 168C 3 2 ) e 3 

— ^— Cs - 1620C 6 ,2 - 3252C 5 C 3 ) £ 5 
+ (98490C 9 - 14598C 5 C4 - 15390CeC3 + 2676C|) £ 6 + 0{e 7 ), 

1 



(l-2e) 



where 



6,2 



E ■ 

ni>7i2>0 
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(85) 



(86) 



(87) 



Finally, eqs. (|85H87[) together with eq. (|72|) lead to the following (independent from our calcu- 
lations) result for MI M 33 which is not only in full agreement to eq. ([FT]) but also includes two 
more terms in e: 

367253 97982^ 11038 , 2 
M 33 ,4 = — + —3— Ca - —3— Cl + 7987C4 - 1966CsC4 

+ 22750C5-3914C3C5 + ^^C6 - 4860C6,2 + 22676C 7 +^^C8, (88) 

2073833 580022 ^ 66370 2 47918 , 
M 33 ,5 = 3 + ^3-^— C3 + -9- C 3 3 + 48991C 4 - 11038C 3 C4 

+ 123766C 5 - 7828C 3 Cs - 35031^5 + ^^Ce 

_ ^ (>Ct _ 9r20(6 , 2 + 1 03 8 38C 7 + H^l C8 + ^ C9 . m 

In the same way we have successfully checked all other simple Mi's. In addition, for all of them 
we get two extra terms of the e-cxpansion. They look similar to the ones listed in (|88I89[) and 
include, in addition to £ 3 — £9, only £6,2- 

5.2.3. Complicated integrals 

We start from diagrams M52 and A/43 which arc relatively simple as they could be expressed 
through the e 1 and e 2 extra terms of the basic three-loop non-planar integral No, namely: 

M 52 = -N (e) = ^ + N QA +N . 2 e + O{e 2 ), (90) 

e e 

M 43 = G (l,2 + 3e)iV (e) = -^ + +25C5-^ + (^-%-75C 5 ) e + O( e 2 ).(91) 

e 4 4 4 



Thus, 

The coefficient TVo.i was known since long from calculations of the five-loop /3-function in the 4> A 



^o,i = -80C 5 + 68C 3 2 + 50C 6 , iV , 2 = -272C 3 2 - 200C 6 + 204C 3 C4 +450C 7 . (92) 



model in [76j. The second coefficient Nq^ was first computed with the GaC method and presented 
in [77J. Its completely independent calculation (through fitting a high- precision numerical result 
with an appropriate analytical ansatz) was performed in [78j |. Needless to say that the results of 
and [3] are in agreement with eq. (|92j) . 

All other complicated integrals (except for convergent integral M44, whose value, ^C7j a t 
D = 4 was also analytically found in [76J) have not been known with sufficient accuracy before 
our calculations. Note, that for a given (four-loop) master integral Mi one needs to know only the 
5 + pi first terms in its e-cxpansion. (For accounting purposes we assume that every expansion 
starts from even if the corresponding term drops out from a specific MI.) Among them first 
5 + p l — 1 (that is all except for the last one) are in a sense easy as they all could be analytically 



found by well-known methods based on Infrared Rearrangement (see, e.g. [19( where the issue 
was spelled out on the example of massive four- loop tadpoles). At any case they all are very well 
checked in the course of the renormalization procedure. 



6. Perspectives 

6.1. Five-loop master integrals 

As we have seen from the discussion in the previous section all complicated and even simple 
four-loop Mi's have been completely expressed in terms of very simple watermelon-like primitive 
p-intcgrals. In fact, the reduction method based on the 1/D expansion of coefficient functions of 
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Mi's is in general applicable for any number of loops in the A-Problem. The GaC symmetry is 
also not limited by number of loops of p-integrals. 

Thus, if the five- loop A-problem is solved, then the five- loop B-problem can also be solved 
in the following sense: the identities stemming from the GaC symmetry will express all five-loop 
Mi's in terms of significantly smaller set of p-integrals. But which exactly set? At present nobody 
knows for sure. But one could certainly expect that: 

• in general the five-loop master p-integrals will contain irrational terms of weight not higher 
than 9; 

• the "small set" of five-loop integrals will include ones primitive as well as those expressible 
in terms of the generalized F-function. 

As both types of the integrals could certainly be analytically evaluated up to the weight 9 



2lL |75| we conclude that the five-loop B-problem should be analytically doable. Moreover, we 
believe that the GaC symmetry + reduction provides the simplest way of analytical solution of 
the five-loop B-problem. 

6.2. General case 

At first sight the applicability scope of the GaC method is rather limited and amounts exclu- 
sively to the massless propagators. Indeed, the heart of the method is the existence of relations 
between integrals of different topologies beyond those provided by the very integration by parts. 
We are not aware about existence of such relations in general case except for a one: finiteness. 

Indeed, two finite integrals are in certainly equal to each other with accuracy O(e ) irrespec- 
tively on their topologies. As a result IBP relations will provide some partial information about 
the values of corresponding master integrals. Unfortunately, the information proves to be rather 
limited^]. 

Indeed, let us consider, as a simple example, eqs. (|24p . Without any use of GaC symmetry one 
could, obviously, write 

7V = io + 0(e°), N Q = Nx + O(e ), N = N 2 + O(e°), N a = O(e ). (93) 

After reduction to masters and solution of the resulting equations we arrive to the same results 
(I40M3|) but with the e-accuracy downgraded by one for every master integral. This is certainly not 
enough to solve the three-loop problem: no new information is obtained for the most complicated 
non-planar master integral Nq. 

The reason for the failure is quite clear: the equations ([Ml) do not provide, in fact, any 
constraints on the value of No as they could be equivalently rewritten as follows: 

L Q = 0(e°), Nx = O(e ), N 2 = O(e ), N a = O(e ). (94) 

Repeating the same exercise at four-loop level we will arrive to a similar conclusion: without 
any use of the GaC symmetry one could find for every master integral Ms all except for the last 
one (that is 5 + p l — 1) of its e-expansion. Again finiteness only, without the GaC-symmetry, is 
not enough to solve the four-loop Problem. 

On the other hand, any L-loop MI multiplied by a one-loop scalar massless propagator is, 
obviously, a (L+l)-loop MI (compare, for, example, T\ and P 2 ). Within the G-scheme framework 



21 Nevertheless, there are cases for which even this limited information is enough, see e.g. I79t, |8' 
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the values of the integrals are trivially related by a factor e. Thus, at two- and three-loops we get 
the following identities 

Ti = eP 2 , T 2 = eP A 

and 

N =eM 52 , Li=eM 32 , Pi = £M 24 , P 2 = eM 23 , P 3 = eM u 

respectively. 

Using "downgraded" eqs. (|58p and (jS(J|) we find the same results for both non-primitive Mi's 
iVo and L\ as in eqs. (|40j) and (|41j) but with the deepness of the e-expansion increased by one. 
Thus we arrive at a truly remarkable conclusion: by merely reducing finite four-loop propagators 
to the master integrals and without any use of GaC symmetry we could not only completely solve 
the three-loop Problem, we even can get one more term in e-expansion of every non-trivial master 
integral! 

It remains to see how predictive is this trick of finding master integrals for cases with other 
patterns of external momenta and masses. But its is absolutely clear that at least some useful 
nontrivial information can be obtained along these lines. 

7. Even zetas 

7.1. Four and five loops 

In addition to six remarkable features of four-loop master p-intcgrals listed in subsection 15.11 
there exist the seventh one, probably most remarkable. Indeed, a scrupulous inspection of eqs. (|47I 
[73])) demonstrate that all their right hand sides do depend on only the following three combina- 
tions of zetas: 

Cs = Cs + yCi - — Ce, Cs = Cs + yCe and ( 7 . (95) 

This simple fact has far reaching consequences. Indeed, a little meditation on (|95p leads to the 
following statemenj^l: 
Theorem 4. 

1. Any finite at e — > p-integral does not contain even zetas {C, 2n \ n > 2} in the limit of e 0. 

2. Any finite at e — > combination of p- integrals like 

j 

with the coefficient functions being functions (not necessarily finite at e — >• 0) with purely 
rational coefficients Cy , will not contain even zetas in the limit of e — 5- (while odd zetas 
{(,2n+i | n > 1} are expected and indeed appear in general). 

3. Let F(e) be any renormalized (and, thus, finite in the limit of e — > 0) combination of any 
p-integrals. The sole source of possible even zetas in F (0) is the appearance of zetas (not 
necessarily even) in the renormalization factors involved in carrying out the renormalization 
of F. 

The third point suggests a clean explanation of an old puzzle of pQCD: the absence of even 
zetas in the Adler function of pQCD at order a^, D^{q 2 ). Indeed, the function at this order is 

(i) a finite combination of four-loop p-intcgrals; 

(ii) the corresponding renormalization is done with the help of charge coupling renormalization 
which docs not depend on any zetas at the order required. 



By any p-integral below we understand any one with number of loops less or equal to four. 
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As a direct consequence of (h) the function D^(q 2 ) should not depend on £4, £g and £3^4 



This is indeed the case [61|, [62]! At first glance, this explanation of this old puzzle generates 
another one: why D^{q 2 ) does not include terms proportional odd zetas with weight larger than 
five (which is expected in general for any combination of four- loop p-integrals)? No, it is not a 
puzzle since long. It is a well-known factl 23 ! that the (L+l) loop Adler function (in massless QCD) 
could be completely expressed through L-loop p-integrals. 

Exactly the same reason explains the absence of even zetas in the four-loop contribution to the 
Bjorken sum rule [io^ . 

The problem of why the five-loop 0(a^) Adler function is also free from even zetas should 
be possible to solve by extending the above reasoning by one loop higher. The "only" missing 
ingredients — a property of five-loop master p-integrals analogous to f|95|) . We hope to come back 
to the subject in future. 

7.2. Three loops 

In fact, Theorem 4 was proven for the three- loop p- integrals in an early work by Broadhurst 
[72I with the help of essentially equivalent (though, to our opinion, somewhat more complicated) 
considerations. This is certainly enough to explain the absence of £4 in the a 2 s contribution to 



Adler function 81 1 and in the a% ones to the deep inelastic sum rules found in 82j, |45| . This 
is because all these quantities are naturally expressed through some combinations of three-loop 
p-integrals with purely rational coefficients. 

Note that, the three- loop version of Theorem 4 (t2| is not enough to explain the absence of 
even zetas from the contribution to Adler function and from the four-loop result for the QCD 



/3-function [83|, |79(. The problem is that, to the best of our knowledge, it is not known whether 
one could find a representation, say, the four-loop contribution to the QCD /3-function in term of 
a finite combination of the three-loop p-integrals with coefficients free from any zetas. The same 
is true for the Adler function. 

On the other hand, we do agree with [7^] that the four-loop QED /3-function should be free from 



£4 (in agreement with explicit calculations of 84[ ) as it can be expressed via a finite combination 



(with purely rational coefficients) of three-loop p-integrals [85[ . 
8. Discussion 

There are various points deserving further discussion in connection with the algorithm of 
evaluation of MPs elaborated in sections H] and [5) 

• The results discussed in the present paper have been indispensable for the long-term project 



of computing the cross section of e + e~ annihilation into hadrons at order a* in QCD 28j, 
While they were first obtained in 2003, their publication had been postponed in favour of 
the faster completion of the main project. 

• Recently, a definite class of massless p-integrals was proven to be expressible in terms of 
the multiple zeta values for all orders of expansions in D — 4, and a direct method of their 



evaluation was suggested ( [86l. l87j). 

In our case (four-loop p-integrals) it predicts that the result is expressible through £ ra up to 
n = 7, as confirmed by our calculations. Unfortunately this method in its present form is 
applicable only for integrals with number of internal lines equal to doubled number of loops 
plus one, so the most complicated four-loop Mi's seem to be unreachable. 



23 First, probably understood on the example of the three-loop C(cr 2 ) Adler function [81 
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The heavy use of identities between Feynman integrals coming, eventually, from IBP relations 
is not a unique feature of our approach to the evaluation of Mi's. It is of interest, that three 
other, quite different and in a sense more general approaches would also be impossible 
without intensive use of the reduction of Feynman integrals to masters. We mean (i) the 
method of differential equa tionJ^l. (ii) the use of difference equations p, |98j and, at last, 
very new me thod [99M1O ll | based on recurrence equations with respect to the space-time 



dimension D 102] 



Finally, we want to mention two popular and powerful methods of evaluation of Mi's which 
do not use directly the IBP reduction. The first approach is based on the Mellin-Barn rep- 
resentation. The ear ly applica tions of Mcllin integrals to evaluation of FI's were performed 



in pioneering works 103l Il04j |. Currently it is an actively developing field, for a review see, 
e.g. Q9- 

The second method — the so-called sector decomposition — was orig inally used as a conve- 



nient theoretical tool for the analysis of convergence of FI's 11051 1106B . First applications of 
sector decomposition for evaluation of F I's w ere considered in 107Hl09j . The current status 
of the method can be found in a review 11101 . 



While the tests of our results described in subsections 15 .2 . 1 1 and 15 .2 . 21 leave no room for doubt 
as for the cases of trivial and simple groups of Mi's, it is not true for the most difficult group 
of complicated integrals: for this family of thirteen Mi's only three had been directly checked 
in an independent way (see subsection 15.2.3)) . It means that if a master integral from the 
remaining ten integrals were assigned a wrong value, it would change in all probability all 
physical results obtained with the use of these Mi's since 2004. 

Let us, therefore, discuss a little bit further the important issue of the correctness of these 
Mi's. The method of computing master p-integrals described in sections [3] and [3] heavily 
uses both the GaC symmetry and the procedure of reduction of four-loop p-integrals. The 
latter is the most complicated part of all the calculation as it requires, first, careful computer 
algebra programming and, second, large-scale calculations. Thus, an independent check of 
the ten remaining most complicated Mi's would also provide us with a quite strong, though 
non-direct, evidence for the correctness of the reduction algorithm we use and its FORM 
implementation. 

Fortunately, s uch an independent check of all complicated integrals has been very recently 



performed in [lllj with the use of the sector decomposition, where not only all results have 
been (numerically) confirmed with better than l%o accuracxr^l but also one extra term in 
these e-expansions have been computed. 



9. Summary and Conclusions 

In this work we have presented an algorithm for the analytical evaluation of all master integrals 
which appear in the process of reduction of massless dimensionally regulated Feynman integrals 
with one external momentum (p-intcgrals). The algorithm is based on the glue-and-cut symmetry 



24 Starting from early works I88h96I the method has developed into quite a powerful technique. For its modern 
status and further references, see the review [97I . 

25 We mean the accuracy for the most complicated last 0(e Pi ) term in comparison with the exact results listed 
in eqs. J50159I) ■ the accuracy of simpler terms of order e l with —4 < i < pi) is significantly higher. The typical 
accuracy of the 0(e Pi + 1 ) term (which is necessary only in evaluation of five-loop master p-integrals) was also about 
l°/oO- One should keep in mind that the latter accuracy is an estimate given by the MC-integrator and as such it 
is not always reliable. 
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[l[ which is an unique and very specific property of such integrals valid irrespectively of their com- 
plexity (number of loops). In addition to the symmetry the algorithm heavily uses the reduction 
procedure. 

It has been demonstrated that the algorithm works flawlessly for the case of the three-loop 
p-intcgrals (successfully reproducing well-known thirty years old results of 21[) and four- loop p- 
intcgrals. In the latter case it produces explicit analytical results for all master integrals, major 
part of which are new. 

Together with Theorem 2 and the 1/D method of reduction of p-integrals 1J, l23[ the algorithm 
guaranties that the UV counterterm of any five-loop diagram can be calculated within the MS- 
schcme in terms of rational numbers. £3, £4, £5, £g an d C7- This implies the analytical calculability 
of the /3-functions and anomalous dimensions of fields and composite operators in an arbitrary 
model at the five-loop level. 
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